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Abstract 

The s-wave correlated 7r-7r strength distribution in cold nuclear matter is in- 
vestigated. The in-medium single pions are renormalized by means of the usual 
p-wave couplings to nucleons and A's. Two groups of models for the tt-tt interac- 
tion are studied. The first one uses phenomenological models which reproduce the 
phase shifts but do not respect the soft-pion theorems. With increasing density 
a strong accumulation of strength in the subhreshold region is observed. Slightly 
above saturation density there even occurs an instability with respect to s-wave tt-tt 
pair condensation. These effects are strongly moderated or disappear altogether if 
constraints from chiral symmetry are built into the tt-tt interaction. 
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1 Introduction 



Correlated two-pion exchange is a crucial piece of nucleon-nucleon attraction at interme- 
diate range 0. On the other hand pion-nucleus phenomenology implies that single-pion 
propagation is strongly modified in a nuclear medium through couplings to nucleons and 
A's [||. It is therefore interesting to ask to what extent pion-pion correlations are changed 
in a nuclear medium. In a first investigation of this kind, Schuck et al. have shown that, 
indeed, the influence of the medium on the it-it correlations may be extremely strong in 
the scalar-isoscalar (cr-meson) channel. Subsequently also the vector-isovector (p-meson) 
channel was investigated showing more moderate effects. The strong influence of the 
medium in the I=J=0 channel arises from the fact that there is rather strong attraction 
between the pions at low energies close to the 2m n threshold. The in-medium suppression 
of the single-pion kinetic energy through coupling to A's and nucleons then leads to a 
threshold accumulation of two-pion strength || or even the occurrence of quasi-bound 
states. For moderate compression it was found in || that the system may become un- 
stable with respect to tt-tt pair condensation. This instability is much softer than that 
for single-pion condensation discussed long ago ||. Such spectacular effects call for more 
exhaustive investigations which are being carried out in the present paper. We first repeat 
the calculation of |5[] using a quite different, but very well tested, phenomeno logical tt-tt 
interaction (the Jiilich model J/J). This model also predicts the instability (with slightly 
lower) critical density and thus confirms the findings of ||. 

A dramatic reshaping of the tt-tt strength distribution inevitably would have strong and 
potentially disturbing effects on the nuclear equation of state since it implies a very long- 
range scalar interaction near the critical density. One has therefore to look for physical 
effects which possibly could oppose the softening of the tt-tt mode. One serious defect of 
phenomenological tt-tt interactions is that they often do not respect constraints from chiral 
symmetry which are known to govern the near threshold behavior ||. This may not be 
crucial for energies well above the 2m n threshold but for subthreshold effects it certainly 
is. We therefore investigate whether the implementation of the soft-pion theorems can 
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hinder an extreme softening of the tt-tt distribution. In a first attempt we use the linear 
sigma model adjusting the parameters such that the 1 = and 1 = 2 scattering lengths 
as well as the low-energy phase shifts are reasonably well reproduced. Indeed it is found 
that this model yields a tt-tt interaction which, when continued off-shell below the 2m,r 
threshold, builds up a strong repulsion. The in-medium calculation of the tt-tt distribution 
then shows very little subthreshold strength and no pair condensation. The linear sigma 
model does not include the p-meson which is a crucial ingredient for the s-wave attraction 
of the Julich model. To study the influence of the p-meson (under the constraints of 
chiral symmetry) we examine, in a second step, the Lagrangian || of Weinberg, which 
incorporates the p-meson as a gauge boson in the non-linear a model. Also this model gives 
no significant subthreshold strength, emphasizing the generic features of chiral symmetry. 
The most crucial point is the correct behavior of the scattering lengths in the chiral limit 

— > 0. We have therefore modified the Julich model to impose the correct chiral limits. 
This leads to the same conclusions suggesting strongly that a tt-tt s-wave pair instablity 
in dense nuclear matter does not exist. 

The paper is organized as follows: To set the stage in section 2 our model for the 
renormalization of the single-pion propagator is outlined. It includes p-wave pion-nucleon 
and pion-A interactions with explicit treatment of the particle-hole continuum and the 
in-medium A width. Section 3 is devoted to the calculation of s-wave correlated pions 
in nuclear matter by using the two phenomenological models of Johnstone-Lee |9j and 
the Julich group J7|. Both predict a pairing instability. In section 4 the constraints from 
chiral symmetry are then included in the framework of the linear sigma model, the gauged 
non-linear sigma model and an extended version of the Julich model, as discussed above. 
Section 5 gives an overall discussion of our findings and concludes with an outlook. 
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2 Single-Pion Propagator 

Let us begin with a brief description of the model for renormalizing the single-pion propa- 
gator in nuclear matter. Including nucleons and A-isobars as the most important degrees 
of freedom and using p-wave dominance in nN scattering the system is described by the 
Hamiltonian 

H = T* + T N + T A + V" f + V? A + V™ + V»Z + V£? , (1) 

where T w , T N and T A denote the kinetic energy operators and V^ N , V^ A are the ttNN and 
7riVA interactions, respectively. Short-range interactions among the nucleons and isobars 
are described by the potentials , V$ A and V AA . The in- medium pion propagator 
can then be extracted from a set of coupled-channel integral equations (Fig. 1) 



D n = V w + V w V wN r Nn + D w V wA r An , 

N , N-i/N n _,_ , N j/N N N , , Ni/NNN 

T Nn — AN V Nn IJ n -T XN V NN T Ntt "+" AN V NA T Aw ) 

T N _ ^NyN D , NyNN N , JVyiVJV AT /q\ 

A-7T — AA ^Att-^vt -I- Xa 17 AA r A7r + AA K A7V r jV7r 5 l Z J 



= - t— ( 3 ) 



where £)° denotes the bare pion propagator 

1 

a; 2 — u;^ + i?7 

with = q 2 + ml, while xjv an d Xa are the particle-hole and A-hole Lindhard functions, 
respectively. Parametrizing the fermionic interaction potentials by Migdal parameters it 
is then straightforward to cast the pion propagator in the following form: 

Dtt = — - (4) 

with the self energy given by 

y ( . n ^ _ 7 2 n pfe + Tl Ah + W h U Ah (2g' NA - g' NN - g' AA ) 

A ' qj q (i + <?Wn^)(i + <? AA n^) - «$ A rM^ ' l ° J 

where 

IPV, q) = -4^|^^ / ^G%(p)G%( P + g) (6) 



and 

q) = - / J^G%( P )Gl(p + g) (7) 

with T n (q) being the usual monopole-type form factor at the 7rNN and 7rNA vertices 
with a cutoff of 1.2GeV/c. Using the nucleon and A Green's functions, G% and G A , 
the integrals in eqs. (2.6) and (2.7) can be easily expressed in terms of the ph and Ah 
Lindhard functions [fj, [TIJ . Due to the high-energy cut in the A-h bubbles, Fermi motion 
effects in eq. (5) can be ignored (which we have verified numerically), and we can put U Ah 
in the simple formf] 

K J 9 ml ^\^- eAh ( q ) + zIr A u + e Ah (ci)i VJ 



with 



eAh(q) = M A -M N + ^-. (9) 



To account for the finite life time of the A-isobar, one has to consider a complex A-energy 
with an in-medium A-decay width, T A . In this study, we shall include both Pauli blocking 
and relativistic corrections as well as a term which accounts for higher-order in-medium 
corrections such that 

r A = r« + ri 2) . (10) 

is taken from reference [IT . For F A we adopt the following phenomenological ansatz: 

ri 2) = 48^ ( 1 + ^L.) A( W ), (11) 



Po V 380 — m v . 

where A (a;) is a cutoff function, which suppresses T A smoothly below and above some 
400 MeV. This choice is motivated by the observation that the inelastic part of the A-self 
energy in the two processes: A — > 7riV and A — > 'jN turns out to be very similar. This 
seems to indicate that T A is insensitive to the nature of emitted particles other than the 

nucleon and one is encouraged to use it for off-shell pions as well. Although somewhat 
t Notice a misprint in the formula for in the reference ||. The results were produced 
with the correct ansatz. 



4 



rough, the ansatz (11) aims to account for the well-known p 2 behavior in the inelastic part 
of the 7r-nucleus optical potential a la Kisslinger [|1^]. It is also important to keep in mind, 
however, that this contribution is not sufficient to reproduce quantitatively the absorptive 
part in the optical potential. One therefore has to consider also contributions from the 2p- 
2h excitations induced by the pion. Following the exhaustive analysis of the Lyon-group 



13, uM, we adopt the following phenomeno logical expression for the imaginary part of 



the 2p-2h pion self energy: 

/mn 2 P- 2/i ( W ,9) = -0.3lt(q)(-^) 2 A( w ) . (12) 

Po 

(2) 

The coefficient 0.3 is fitted such that the combined effects of T A in the RPA approximation 
and the 2p-2h contribution properly accounts for the empirical absorptive part of the 
7r-nucleus optical potential at threshold. The real part of the 2p-2h pion self energy, 
although straightforward to evaluate via dispersion relations, will not be considered since 
its influence is minor in the present study. For practical purposes we absorb n 2p_2/l into 
the A-hole polarization propagator as 

This concludes the outline of our in-medium renormalization of the pion propagator. 

3 7T-7T Correlations in Nuclear Matter: Phenomeno- 
logical Models 

3.1 The Johnstone-Lee Model 

To study the tttt correlations in dense nuclear matter, we first introduce a simple phe- 
nomenological model for the scattering of two pions in vacuum. This model is based on 
the separable potentials of Johnstone and Lee |9j . For two-pion CMS energies below ~ 1 
GeV these give a realistic description of the tttt phase-shifts in low partial waves. Aside 
from numerical advantages, the separability renders simple and transparent expressions 



for the cr propagator. In the next section we shall use the more sophisticated meson- 
exchange model by the Jiilich group, which is able to extend the description above the 
KK threshold. 

In the model of Johnstone and Lee, the interaction kernel for tttt scattering in the 
J p = + , I = channel is of the form 

V^iq, q', E) = 2n 2 gl2u q 2uj q ,v <T {q)Dl{E)v (7 {q')i (14) 

where 



Ha 

denotes the finite-size otttt vertex and 

= (16) 

the bare cr-propagator. The coupling constant g a , the cutoff q a and the bare a-mass 
M a are free parameters, adjusted to the experimental 1=0-, s-wave phase shift and the 
scattering length. The T-matrix in "back to back" kinematics then takes a very simple 
form 

T™(q, q>, E) = 2<K 2 gl2u q 2uj q ,v a ((l)DAE)v a (q'), (17) 
which is the same as (14) except that is replaced by the dressed propagator 

= g-M? + ' <18) 
Here S CT is the u-meson self energy due to tttt rescattering processes involving s-channel 
pole graphs. Explicitly 

- %L a (E) = 2n*gl J ^^{q)^ ^(E -u,q). (19) 

In nuclear matter, modifications of the T-matrix enter through the dressing of the pion 
with ph and Ah excitations (Fig. 2). This amounts to replacing D® with the renormalized 
propagator D n given in sect. 2. In previous work M, we have used a simple two-level model 
by coupling the pion only to the A-hole branch. As a result, a strong accumulation of 
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strength near the 2m n threshold was found in the sigma-mass distribution. For densities 
somewhat above saturation density even a two-pion bound state emerges. The shift of 
the a strength to low energy can be easily understood as being due to the softening of the 
pionic branch induced by the well-known level-mixing mechanism. In the present paper 
we can go beyond this simplified treatment by incorporating the particle-hole branch, 
including Fermi motion. It is clear that any strength accumulation below the two-pion 
threshold will depend sensitively on the energy distribution of the ph excitations, thus 
Fermi motion and Pauli blocking effects are crucial. Secondly we include a realistic treat- 
ment of the A-decay width as well as absorptive effects as described in sect. 2. With 
these extensions the energy integration in eq.(19) can no longer be done analytically as it 
is the case in the pure two-level model of [[|. Obviously a numerical integration over the 
entire range of energy is not a trivial task. To overcome this difficulty we first compute 
the imaginary part of S CT given by: 

ImE a (E) =g 2 a J dqqHu 2 q vl{q)ImG^{E, q) (20) 

with being the two-pion propagator (the convolution of two single-pion propagators). 



Its imaginary part can be computed by using the following expression [g, Ts 

1 r E 

ImG nn (E,q) = / dLjImD w (uj,q)ImD n (E — u,q) , (21) 

7T Jo 

which simplifies the numerical integration considerably since off-shell integrations can be 
avoided. The real part of S CT is computed by means of the usual dispersion relation, that 
is 

where V denotes the Cauchy principal value. To obtain results we have to specify the 
Fermi liquid parameters, which represent the short-range repulsion in the nucleon and A 
interaction potentials as discussed in sect. 2. We have chosen the following values: 

9' NN = 0.8, (23) 
9na = 9'aa = °- 5 



based on empirical knowledge from Gamow- Teller systematics \T6\ as well as information 



from pion double-charge exchange reactions [IT). In Fig. 3 we contrast results from the 
simple two-level model with those including isobar-decay and absorption effects. As ex- 
pected, these modifications lead to a broadening of the bound state but have a minor 
effect on the overall distribution. The situation changes if the particle-hole branch is 
included. In this case the two-pion state couples to the particle-hole continuum, resulting 
in an enhancement of the decay width below threshold as can be seen from the upper 
panel in Fig. 4, where the imaginary part of S CT is displayed. The width strongly increases 
with density.In the real part of D a (lower panel in Fig. 4) one observes a rapid decrease of 
the 2ir pole towards zero energy, signalling the pairing instability. The critical density is 
slightly above 1.3 po- This is significantly lower than the threshold for pion condensation, 
which in our model occurs at 3.6po- The reason is the extra attraction provided by the 
s-wave 7T7T interaction. The approach to pair condensation is also seen clearly in the a- 
strength function (Fig. 5) where we notice that, with increasing density, a resonance below 
threshold develops which becomes narrower and absorbs practically the entire strength 
close to the instability. 

3.2 The Jiilich Model 

In order to check to what extent the results depend on the specific form of the irn inter- 
action we have repeated the calculations for the model of Lohse et al. 0. Here V n7T is 
constructed from meson-exchange including the most important s- and t-channel contri- 
butions. Besides the a pole term, the t-channel p exchange is very important for s-waves 
(Fig. 6). In addition, the model couples the KK channel, which becomes crucial for 
energies > 1 GeV. To derive an equation for the T-matrix, a Blankenbecler-Sugar (BbS) 



reduction [18| of the 4-dimensional Bethe-Salpeter equation is employed maintaining co- 



variance The T-matrix for given angular momentum J and isospin I is then obtained 
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as 

2 poo 

Tg(E; q, q') = V£{E; q, q') + + £ / rffc A; 2 4 W 2 V; J /(£; g, fc) G° X (E, k) Tg(E; k, q') , 

A=l J ° 

(24) 

where the indices /i, A=l,2 label the two different channels hit or KK, respectively (also 
io\ = m\ + k 2 with m\ = m^^rriK)- In the BbS form the vacuum 7T7r-/i^^-propagator is 
given by 

G\(E, k) = — F2 ] 2 . (25) 

uj k E 2 -Auol + iT] 

This model gives a good description of the phase shifts and inelasticities up to ~ 1.5 GeV. 

For the present problem it is most relevant that the Jiilich model has a different off-shell 
behavior than the Johnstone-Lee model which could potentially alter the subthreshold 
behavior of D a since the pions are far off shell. 

Because of the i-channel interaction the calculation of the a propagator is more in- 
volved than in the separable model. The way we proceed is as follows: in a first step we 
calculate the full off-shell T-matrix with the two pions embedded in nuclear matter. In 
analogy to sect. 3, we make use of a dispersion relation, i.e. we calculate ImG^n using 
eq. (21) and evaluate the real part as 

ReGM^-VrdE-l^^l. (26) 

The T-matrix equation is then solved numerically with the replacement G® n — > G^ (the 
KK propagator remains unmodified) employing the Haftel-Tabakin method [^(J . The in- 
medium T-matrix, T°°, then serves to dress the bare cx-propagator given by eq. (16) with 
a bare mass of M a = 2.2 GeV (in free space the renormalization of D®(E), generated by 
the iteration of V, leads to the /o(1400)-resonance). The medium-modified cx-propagator 
then takes the following form 

2 poo 

D a {E) = D° a (E)+D°(E) £ / dqqHuy^E^q) G x (E,q) v X a (E,q) 

x=i J0 

2 poo 

+ £ / dqq 2 4u 2 vZ(E,q) G,{E,q) 



x £ / dq'q' 2 4:UJq,T™(E; q, q>) G U (E, q>) <(£, </)] , 

(27) 

where v^(E,q) denote the vertex functions for the s-channel transitions tttt, KK — > a. 
The second term on the right hand side of eq. (28) accounts for the single im-/KK- 
bubble, whereas the third term sums all possible tttt- and fT.ft'-rescattering processes. We 
should emphasize that it is necessary to include the contributions from the KK channel 
which, off-shell, enters below the KK threshold at 2m k- Neglecting this channel results 
in an unphysical pole at the bare mass M a . For energies E > 2m k we have introduced 
a small width for the kaons to avoid numerical instabilities (moderate changes in the 
width parameter affect our results for D a (E > 2mx) only slightly; this region is of minor 
interest in the present discussion anyway). Fig. 7 displays the results for density values 
p/Po — 0.5,1.0 and 1.2. They are very similar to those obtained with the separable 
model (Fig. 5). The critical density for pairing is somewhat lower at p/po ~ 1.3, which 
can be traced back to the somewhat larger s-wave attraction in the Jiilich model. This 
result seems to suggest that any tttt model, which reproduces the phase shifts, leads to an 
instability at densities slightly above saturation density. 



4 tt-tt Correlations in Nuclear Matter: Chiral Models 

As mentioned in the introduction, the models used in the previous section do not respect 
the constraints from chiral symmetry, i.e. the soft pion theorems which imply a linear 
decrease of the s-wave tt-tt 1 = and 1 = 2 scattering lengths as m^ — > 0. We will now 
elaborate in more detail on these properties and then apply various chiral models below. 

It is possible to relate the I=J=0 tt-tt scattering amplitude at the soft point, defined 
through the Mandelstam variables s = m^,t — 0, u — m%, to the tt-tt sigma commutator 

M 00 (s = mlt = 0,u = ml) = ^, (28) 
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with Etttt given by 

=< 7r*(k = 0) | [Q 5 , [Qs, #]] K(k = 0) > eomeeted , (29) 

where Q5 is the axial charge and H the full strong-interaction Hamiltonian. The (discrete) 
normalization is such that 

< 7r i (k)|7r i (k') >= 2cj k <5 kk /. (30) 

The derivation of this result only requires the standard reduction formula and PCAC 
and follows exactly the same lines as the derivation of the 7r-nucleon sigma term (see for 



instance standard text books such as PHI ). An explicit calculation shows that the double 
commutator appearing in eq. (29) is just the piece of the Hamiltonian which explicitly 
breaks SU(2) x SU(2) chiral symmetry. Hence, in an effective theory with pions, this 
yields 

= J d^-m\ < vr(k = 0)| : $?(x,t) : |vr(k = 0) >= m\. (31) 

where is the canonical pion field with isospin i. We thus obtain a low-energy theorem 

2 

TYl 

M 00 (s = mlt = 0,u = ml) = -f. (32) 

The reader may notice that this result is valid up to m\ corrections as it is implicit 
in PCAC. Eq. (32) should be contrasted with the Weinberg result || at the physical 
threshold which gives 

M 00 (s = Ami, t = 0, u = 0) = (33) 

J IT 

and indicates a sign change when going off-shell. However, the soft-point amplitude is 
not directly relevant to the tt-tt scattering amplitude where the Mandelstam variables are 
to be chosen as 

s = E 2 , 

t = 2ml ~ Z^qUqi + 2qq', 

u = 2ml ~ 2uj q uj q i — 2qq' (34) 
11 



with q and q' being the relative momenta of the in- and outgoing pion in the CM frame. 
In the medium it will be possible to reach kinematical conditions which are relatively close 
to the point where the soft pion result applies. For instance at q = q' =0 and E = m n one 
has s = ml, t = u = 0. Hence, one may suspect that the soft-pion constraint will generate 
repulsion below the threshold which is completely absent in the phenomenological models 
of Johnstone-Lee and the Jiilich group. 

4.1 The Linear Sigma Model 

The linear a model does obey the soft-pion theorem eq. (32). At the tree level, the 1=0 
amplitude is given by the Born amplitude (m a is the mass of the a meson) 

Mr(M,«) = ^^(3^4 + ^4 + ^4) (35) 
j£ s — m~ t — m z a u — m~ 

which satisfies 

M 00 (, = ml,t = 0,u = ml) = (1 - % ( 36 ) 

fi mj 

We also see that, for s = ml, u = t = 0, which can be reached in subthreshold n-n 
scattering, the amplitude is even more repulsive, namely 

ft 1 ml 

We will now use the Born amplitude (35) of the linear a model to construct a seperable 
7r-7r potential. With the kinematical choice (34) for the Mandelstam variables, neglecting 
the t- and u-dependence in the denominators (t,u <C m a ) and adding phenomenological 
form factors to regularize the momentum integrals, we have 

V(E, q ,fl = ^^(3f^4 + 4 ^: 2ml )vMvW . (38) 

The vacuum s-wave mi amplitude is obtained by solving the Lippmann-Schwinger equa- 
tion 

M™(E a a') - V™(F a a) + f ^ Y^llSl^l^ m) 
M (E,q,q)-V (E,q,q)+J ^ 2wfc(E 2 _ 4^2 + ^ ( 39 ) 



M 00 ( S = ml,t = 0,u = 0) = 2^(1 - % (37) 
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With 

m a = l.OGeV, v(q) = (1 + J^Y 1 (40) 

the scattering length is a 00 = 0.23m" 1 and we obtain a good fit to the phase shift up to 
E w 900 MeV (Fig. 8). 

As before, the in-medium scattering matrix is evaluated by modifying the 2-n propagator 
according to eqs. (21) and (26). Rather than displaying the a propagator, as in sect. 3 we 
here show the invavariant amplitude M 00 as in (39) because for some of the chiral models 
a a propagator cannot be defined. However, for the non-chirally symmetric models of 
chapter 3 the imaginary part of the M 00 amplitude has the same shape as lmD a as 
displayed in Figs. 5 and 7. As can be seen from Fig. 9, the density-dependence of the 
strength distribution (upper panel) is rather weak. On the other hand, the real part (lower 
panel) changes sign at low energy reflecting the repulsive character of the tttt potential 
(38) close to zero energy. This behavior prevents the pairing instability. 



4.2 The Weinberg Lagrangian 

To study the influence of the p-meson we now proceed to the non-linear realization of 
chiral symmetry as proposed by Weinberg ||. Before going to nuclear matter, let us first 
briefly recall the model. In the pure mesonic sector, the Lagrangian can be cast in the 
following form 

4 2(1 + ^)2 2 M + iJ 2 ^ ml 1 + ^ > 

Along with the pions, the p-meson is introduced via gauging, maintaining chiral symmetry. 
Here, p M , F^and m p denote the canonical field, the field strength tensor and the p-mass, 



respectively. The coupling constant g p7T7T is fixed by the KSFR relation p2 



2g%Jl = vr? p . (42) 

To lowest order in the expansion of the pion field one can derive the following Born 
amplitude from eq. (41) 

M| =0 = M( =0 + M 2 /=0 + M 3 /=0 , (43) 
13 



with 



Mt\s,t,u) = -l[3(s-ml) + (t-ml) + (u-ml)} (44) 

J 7T 



2 



Mt^s^u) = [(s- u ) + (s-t)} (45) 
ft 



,^r-n/ \ 2m 2 s — u s — t^ . , 

f 2 t-m 2 p u-m 2 p 



In the above expressions we have used the usual metric in which g 00 = l,g* J — — 1. The 
amplitude M[ =0 arises from the contact terms induced by the covariant derivative and the 
symmetry breaking term of the Lagrangian (2nd and 3rd term on the r.h.s. of eq. (41)). 
The other two amplitudes follow from the p-meson exchange and the contact term arising 
naturally from the gauging procedure (last term on the r.h.s. of eq. (41)). As pointed 
out by Weinberg, M^ =0 and M| =0 cancel at the physical threshold such that only the 
amplitude M( =0 contributes to the scattering length. One obtains a 00 = O.lTm^ 1 as 
predicted by Current-Algebra |23| . 

To go beyond the tree level we supplement, as for the linear sigma model, the Born 
amplitude (43) with phenomenological form factors of dipole type, 



F(q) = [ A2 + 4m 'f. (47) 

In contrast to the previous section, where the on-mass-shell identification of the kinematic 
variables has been employed, we choose the Blankenbecler-Sugar prescription |18| (also 
used in Julich model) to reduce the 4-dimensional Bethe-Salpeter equation. In this case 
the two pions in the intermediate state are put on-energy shell. It implies the following 
identification of the Mandelstam variables 

s = (qi + q2) 2 = E 2 , 

t = (qi-q 3 ) 2 = -{q~q'f , 

u = {qi-q±) 2 = -{q + q'Y ■ (48) 

which then determine the energy- momentum dependence of the pseudopotentials V(E, q, q') 
suitable for iteration. 
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With form factor cutoff values Ai « 800 MeV, A 2 « 1300 MeV and A 3 « 2000 MeV 
we can approximately reproduce the experimental s-wave scattering phase shifts up to 
E « 900Me\/. 

As can be easily verified, the pseudopotential (38) of the linear sigma model preserves 
the correct dependence of the scattering length. With the kinematical choice (48) 
this is no longer the case. This is partly due to the fact that the pseudopotential is no 
longer separable. The scattering lengths a 00 and a 02 actually diverge in the chiral limit. 
Furthermore, the repulsion in the subthreshold region, which is present at the tree level, 
does not survive the iteration. In medium this implies again an s-wave instability for 
slightly compressed nuclear matter. 

Thus, inspite of the special care given to chiral constraints at the tree level, one must 
be cautious when iterating. There is, however, a way of preserving the chiral limit to all 
orders. One employs a subtracted dispersion relation to extract the real part of G nn while 
preserving unitarity. We start from the spectral representation of the 2-pion propagator 

G„(s,k) = -lfdJ ^ G *y . (49 ) 

71 J S — S' + 17] 

In vacuum, the imaginary part reads 

ImG nw (s, k) = -— 5(s - Auj 2 k ) . (50) 

In principle, the real part is only fixed up to a real constant. We make use of this freedom 
by subtracting the value at s=0 

ReG nn (s,k) = ReG n7V (s, k) — ReG nn (0, k) 



TT J 



, ImG wn (s,k) s 

ds 51 

s — s' s' 



GUs,k) = -—— r —— 2 . (52) 



The resulting, once subtracted, 2-pion propagator takes the form 

J_ 1 s_ 

uj k s — AujI + it] Auj'I 

which differs from the non-subtracted one (compare eq. (25)) by a factor s/4u; 2 . Notice 
that the imaginary part remaines unchanged, i.e. 

ImG K . K (s,k) = ImG nn (s, k) , (53) 
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thus maintaining the unitarity. Recalling the BbS T-matrix equation (24) with the once- 
subtracted propagator 

T£{E; q, q') = V^(E; q, q') + jf dk k 2 4u 2 V£{E; q, k) G° nn (E, k) — 2 7%{E; k, q'), 

(54) 

one realizes that the subtraction factor s/4u;| cures the divergence of the scattering length 
since, in the soft pion limit, s oc m 2 . By adjusting the form factors (Ai = 1.5GeV, A2 = 
3GeV, A 3 = 4GeV) we obtain a 00 = 0.20m" 1 and a 02 = — 0.031m" 1 which are in perfect 
agreement with the most recent experimental determinations [2~I|. As a result of the 
correct chiral behaviour of the scattering length (which we also checked numerically) the 
'subtracted' M-amplitude again shows repulsion close to E=0 (full line in the lower part 
of Fig. 10). Above threshold, our fit underestimates the phase shifts (note the different 
scale as compared to Fig. 9). This latter point will be improved in the next section. 

To end this section, let us turn to the in-medium M-amplitude calculated by means 
of the once subtracted 2-pion propagator eq. (51). In section 3 it was shown that tttt in- 
teractions without chiral constraints may lead to condensation. In the linear sigma model 
this was not the case as a consequence of a repulsion in the M-amplitude below the 
threshold. The same mechanism is at work when applying a once subtracted dispersion 
relation to the 2-pion propagator of the scattering equation in the BbS framework. As 
one can see from the upper part of Fig. 10, the repulsion in the iterated M-amplitude 
prevents the instability. 



4.3 The Jiilich model with chiral constraints 

To be more realistic in the description of the phase shifts over a wide range of energies, 
we here present a preliminary version of the Jiilich model including constraints from 
chiral symmetry. As described in the previous section, we impose the correct chiral limit 
for the scattering lengths by applying the subtraction method (51) to both the 27r and 
the KK propagator. It turns out that the l/4u;| dependence of the subtraction factor in 
(54) leads to an appreciable suppression of the pseudopotentials V JI in the integrand at 
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small momenta, especially for the contribution from t-channel p exchange. This causes 
a considerable reduction of the attraction in the (00) channel above threshold. We have 
therefore added the 4-pion contact interactions generated by the Weinberg Lagrangian as 
well as an /2(1270) t-channel exchange in the tttt and in the KK system. The importance 



of the latter has already been pointed out in ref.p5|. The resulting fit to the phase shift 
is not quite as good as in the original Julich model, but it qualitatively reproduces the 
important features in all partial waves up to J=2 and E »s l.lGeV (see Fig. 11 for the 
case of the JI=00 channel, which is the one of interest in the present discussion). The 
resulting in-medium M 00 amplitude reconfirms the findings of the two previous sections. 
Even at very high densities there is only a very moderate accumulation of strength slightly 
below 2m 7r , i.e. the system is far from condensation (see Fig. 12). The reason is the same 
as before namely the real part of the nir interaction changes sign and becomes repulsive at 
energies below the 2m n threshold. When the chiral constraints are imposed this feature 
is preserved at the T-matrix level. Nevertheless, we observe a resonance structure below 
threshold (insert in the upper panel of Fig. 12), which corresponds to a quasi bound 2-pion 
state. 



5 Summary and conclusion 

Based on a conventional description of single-pion propagation in the nuclear medium 
and various models of the vacuum nir interaction we have studied s-wave tttt correlations 
in symmetric nuclear matter. When using phenomenological interactions which give a 
realistic description of the phase shifts in low partial waves, we find 7r + 7r~ s-wave pair 
condensation with a critical density only sightly above saturation density, thus confirming 
earlier findings || . The reason is a strong coupling to the ph continuum combined with the 
attractive s-wave tnr interaction. Such an instability would have dramatic consequences 
for the nuclear equation of state. 

One problem with phenomenological two-pion interactions is that they do not obey 
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the constraints imposed by chiral symmetry, most noteably the correct behavior of the 
scattering lengths in the chiral limit m n — > 0. When constructing scattering potentials 
from chirally symmetric Lagrangians (the linear a model and the Weinberg model) which 
extrapolate correctly to zero pion mass, we have shown that the instability disappears 
and only minor modifications are obtained in the subthreshold region. These findings are 
generic and do not depend on the details of the interaction. To be more realistic, we have 
also studied an improved Jiilich model in which the correct chiral limit is imposed by a 
once-subtracted dispersion relation and which gives a satisfactory description of the phase 
shifts up to 1 GeV. Also in this case instabilities are avoided. 

In conclusion we want to argue, that the in-medium s-wave nn correlations crucially 
depend on chiral symmetry constraints, especially the correct chiral limit for the scattering 
length. Once these are implemented in the theory, the density effects are rather small. 
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Figure Captions 



Figure 1: Particle-hole and A-hole bubbles. 

Figure 2: Dispersion equations for the most important pionic branches in nuclear matter at 
p = pq. The solid curves designate the particle-hole continuum. Dashed curves 
are for the 'interacting' A-hole and pion branches. Here the A-hole and the pionic 
branch are decoupled from the p-h continuum. For the in-medium nucleon mass we 
take 80% of the free one. 

Figure 3: a - strength function for p = Pq. Dashed curve is for the pure two level model of Q. 

The bound state is situated, for this density, at E niT = 1.92m n . The solid curve is for 
a two level model with the corrections mentioned in the text, namely the A-decay 
width and the 2p-2h pion self energy. 

Figure 4: (a) Imaginary part of the a-self energy in the pure two level model (dashed curves 
for 0.5p and p ) an d in the complete model (solid curves for 0.5p , Po an d 1.3p . 

(b) Real part of the inverse cx-propagator, D^: 1 , in the full model. The solid curve 
is for 0.5po, the long-dashed curve for p and the short-dashed curve for 1.3po- As 
explained in the text, the quasi pole appears for p = p at = lSm* and then is 
pushed to E — Q.lm^ at p = 1.3p . 

Figure 5: Evolution of the a-strength function in the full model calculated with the separable 
mi interaction. At 0.5po and po, the imaginary part begins and then invades the 
subthreshold region. For 1.3po, it is constituted in a resonant like structure around 
the quasi pole. The latter condenses for p > 1.3po. 

Figure 6: Lippmann-Schwinger equation for the tttt T-matrix in the meson exchange model 
of the Jiilich group. Below E ~ lGeV the dominant contribution is p t-channel 
exchange. 
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Figure 7: a-strength function for the Jiilich tttt interaction using the full model for the one- 
pion propagator. Here the tttt bound state condenses at slightly lower density than 
in the separable model, namely at p ~ 1.3po- The zero near 7m n is caused by the 
opening of the KK channel. 

Figure 8: Fit to the tttt phase shifts in the JI=00 channel using the linear a model. 

Figure 9: In-Medium results for the imaginary part of the invariant M 00 amplitude calculated 
in the linear a model. Solid curve: vacuum result, long-dashed: p = 0.5po, short- 
dashed: p = po, dotted: p = 2p . 

Figure 10: Invariant M 00 amplitude in the gauged nonlinear a model (Weinberg Lagrangian) 
employing the subtraction scheme: full line: vacuum results, long-dashed line: p = 
0.5po, short-dashed line: p = po, dotted line: p = 2p - 

Figure 11: Phase shift fit within the extended Jiilich model 

Figure 12: In- medium M 00 amplitude calculated with the extended Jiilich model: 

full line: vacuum results, long-dashed line: p = 0.5po, short-dashed line: p = po, 
dotted line: p = 2p . 
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